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We consider the expectation value of a chromo- magnetic flux loop in the immediate forward light 
cone of collisions of heavy nuclei at high energies. Such collisions are characterized by a non-linear 
scale Qs where color fields become strong. We find that loops of area greater than ^ 2/Ql exhibit 
area law behavior, which determines the scale of elementary flux excitations ("vortices"). We also 
estimate the magnetic string tension, gm — 0.12(35. By the time t ^ 1/Qs even small loops satisfy 
an area law while substantial magnetic fields remain only at large distances. We describe corrections 
to the propagator of semi-hard gluons at very early times in the background of fluctuating magnetic 
fields. 



Collisions of heavy ions at high energies provide oppor- 
tunity to study non-linear dynamics of strong QCD color 
fields [l] . The field of a very dense system of color charges 
at rapidities far from the source is determined by the clas- 
sical Yang-Mills equations with a recoilless current along 
the light cone [2] . It consists of gluons characterized by a 
transverse momentum pt on the order of the density of 
valence charges per unit transverse area [3 ; this satu- 
ration momentum scale separates the regime of non-linear 
color field interactions at pr ^ Qs or distances r ^ 1/Qs 
from the perturbative regime at pr ^ Qs- Near the cen- 
ter of a large nucleus this scale is expected to exceed 
-1.5 GeV at BNL-RHIC or CERN-LHC collider en- 
ergies, for a probe in the adjoint representation of the 
color gauge group. The classical field solution provides 
the leading contribution to an expansion in terms of the 
coupling and of the inverse saturation momentum. 

The soft field produced in a collision of two nuclei is 
then a solution of the Yang-Mills equations satisfying ap- 
propriate matching conditions on the light cone [4] . Most 
interestingly, right after the impact strong longitudinal 
chromo-magnetic fields Bz ~ 1/g develop due to the fact 
that the individual projectile and target fields do not 
commute O |6] . They fluctuate randomly according to 
the random local color charge densities of the valence 
sources. In this Letter we show that magnetic loops Wm 
exhibit area law behavior, and we compute the magnetic 
string tension. We also determine Wm for finite times 
after the collision from a numerical solution of the Yang- 
Mills equations. Finally, we sketch how the background 
of magnetic fields affects propagation of semi-hard parti- 
cles with transverse momenta somewhat above Qg. 

We consider a spatial Wilson loop with radius R in the 



plane transverse to the beams, 

M{R) = Vexp (^ajde ^A^j 

WMiR) = ^(trM(i?)) , (1) 

where x = i?(cos ^, sin ^), and path ordering is with re- 
spect to the angle 0; in numerical lattice simulations it 
is more convenient to employ a square loop. We com- 
pare, also, to the expectation value of the Z{Nc) part of 
the loop; for a magnetic field configuration correspond- 
ing simply to a superposition of independent vortices the 
loop should equal exp(27rz n/A/'c), with n the total vor- 
tex charge piercing the loop. Thus, for two colors we 
compute 

W^^'\R) = {sgntvM{R)) (2) 

where sgn() denotes the sign function. Comparing ([T]) 
to Q tests the interpretation that the drop-off of Wm{R) 
is due to Z{Nc) vortices, without requiring gauge fixing 
of the SU(A^c) links [7 . 

The field in the forward light cone immediately after 
a collision [4 , at proper time r = y/t'^ — +0, is 

given by A* = ol\ ^ ol\. In turn, before the collision 
the individual fields of projectile and target are 2d pure 
gauges, 

C^m='-^n^^'^l , S^<=^p^, (3) 

where m = 1, 2 labels projectile and target, respectively, 
and IJrri are SU(N) gauge fields. Eqs. ^ can be solved 
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FIG. 1: Expectation value of the magnetic flux loop right after 
a collision of two nuclei (time r = +0) as a function of its area 
A' = AQ^,. We define = (CF/27r) //x^ Symbols show 
numerical results for SU{2) Yang-Mills on a 4096^ lattice; the 
lattice spacing is set by q^ijll — 0.0661. The lines represent 
fits over the range A> A' >2. 
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FIG. 2: Time evolution of the magnetic flux loop after a colli- 
sion of two nuclei (4096^ lattice, ijll — 0.05). From bottom 
to top, the curves correspond to time r x ji = 0, 1, 2, 3, 
where g^ fi ^ 3Qs so that r = 3/{g^/j.) corresponds to about 
r ^ 1/Qs in physical units. 



either analytically in an expansion in the field strength [4] 
or numerically on a lattice [8 . 

The large-x valence charge density p is a random 
variable [12]. For a large nucleus, the effective action de- 
scribing color charge fluctuations is quadratic, 



P"(x)p^(x) 

2/i2 



(4) 



with /i^ proportional to the thickness of a given nu- 
cleus [2]. The variance of color charge fluctuations de- 
termines the saturation scale ~ 5'^/^^ [3]. The coarse- 
grained effective action Q applies to (transverse) area 
elements containing a large number of large-x "valence" 
charges, AA± /i^ - AA^ Ql/d^ > 1- 

In fig. [l] we show numerical results for Wm immedi- 
ately after a collision. It exhibits area law behavior for 
loops larger than A > 2/(5^. The corresponding "mag- 
netic string tension" is cfm/Q1 = 0.12(1). The area law 
indicates uncorrelated magnetic flux fluctuations through 
the Wilson loop and that the area of magnetic vor- 
tices is rather small, their radius being on the order of 
^vtx ^ 0.8/ Qs- We do not observe a breakdown of the 
area law up to A ~ 4/(5s, implying that vortex correla- 
tions are small at such distance scales. Also, restricting 
to the Z(2) part reduces the magnetic flux through small 
loops but (Jm is comparable to the full SU(2) result, if 
somewhat smaller. 

The numerically small vortex size that we find is para- 
metrically consistent with the classical Gaussian approx- 
imation at weak coupling which, as already mentioned 
above, applies for areas AA_\_ g^/Ql- Corrections to 
5'eff of higher order in p [9 as well as due to quantum 
fluctuations [10 of the fields should be investigated in 
the future. 



Since the field of a single nucleus is a pure gauge it 
follows that W^^^ = 1. Thus, if Wm is computed pertur- 
batively one has to subtract disconnected contributions. 
We estimate the connected contribution to Wm after a 
collision of two nuclei by expanding Wm to fourth order: 

WM{a\ + ai,) - W'^^\al) W^^^K) = 

This result does not require an IR cutoff. There is no 
contribution at order ^ A, hence no area law behav- 
ior and no vortices, even though there is, of course, a 
non-zero longitudinal magnetic field even at the "naive" 
perturbative level: 



2 ^2 

(tr5,(r)5,(rO)=4- ^ 



N? 



1 



A 



(6) 

The failure in eq. ([5| to reproduce the area law be- 
havior of Wm is related to the perturbative expansion of 
the Wilson loop and the interchanged integrations over 
the loop and over the Gaussian charge density fluctua- 
tions Q. In contrast, for the exact lattice result we com- 
pute Wm for every individual configuration of sources 
Pi2{^-l) arid average it a posteriori. Any such particu- 
lar configuration contains higher (magnetic) multipoles 
which are eliminated by a "naive" Gaussian approxima- 
tion; this has been discussed in detail in the recent lit- 
erature [11] which has focused on higher multipoles of 
electric Wilson lines, however. The area squared behav- 
ior predicted by the naive Gaussian approximation ( 5j6 ) 
might be visible in numerical measurements of much 
larger loops where the dipole moment should be domi- 
nant. 
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In fig. |2] we show the time evolution of the magnetic 
flux loop after a collision. The magnetic fleld strength 
decreases due to longitudinal expansion and so Wm ap- 
proaches unity. On the other hand, the onset of area 
law behavior is pushed to smaller loops, implying that 
the size of elementary flux excitations or "vortices" de- 
creases; by the time r ~ l/Qs area law behavior is sat- 
isfled even for rather small loops. Since long wavelength 
magnetic flelds remain even at times ~ l/Qs^ it will be 
important in the future to understand the transition of 
Wm to behavior expected in thermal QCD. 

We have also investigated the dependence of the mag- 
netic flux loop in the adjoint representation on its area. 



Wf^^ = = 



-(|trMp-l) , 



(7) 



and found behavior similar to flg. [2| The adjoint mag- 
netic string tension is about two times larger, as expected 
from (jr]). 

The third color component of the longitudinal mag- 
netic fleld is shown in flg. [3j Domain-like structures 
where the magnetic fleld is either positive or negative are 
clearly visible; they lead to the above-mentioned area law 
of the Wilson loop. Also, one can see that in time the 
magnetic flelds become weaker and smoother. 

The magnetic flelds modify the propagation of semi- 
hard modes with pt not too far above Qg. Quantum me- 
chanically, the transition amplitude from a state 
to \xf^tf) is given by a Feynman sum over paths, 

J ds j Vx^ /exp i j dr {mx^ + gA^x^) \ ^ 

J ds J Vx^ ex.p ^ J drmx^ exp(— ctm^) ,(8) 

where x^{r) is a parametrization of the path with the 
given boundary conditions and length 5; and x^ = 
dx^/dr. Here, the area A is that enclosed by a quan- 
tum mechanical path from the initial to the flnal point 
returning to via the classical path; see flg. |4] The 
classical path is obtained by extremizing the action but 
a single path is a set of measure zero. Semi-classical 
paths can dominate the integral only if there is construc- 
tive interference among neighboring paths from within a 
de Broglie distance. On the other hand, destructive in- 
terference of such paths leads to Anderson localization of 
the wave function. 

Hence, up to a numerical factor, the area in eq. (|8| 
should be given hy A ^ s/pt- Integrating over the 
Schwinger parameter then leads to the propagator 



(9) 
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FIG. 3: Color-3 component of the magnetic field F^y{x,y) in 
the transverse plane at time t = +0 (top) and l/g^fJ, ~ l/SQs 
(bottom) for a single configuration of color charge sources p. 



where ctm = 0.12 Ql from above and m is the mass (time- 
like virtuality) of the particle. This expression accounts 
for corrections to free propagation and could be useful 
for studies of the dynamics of the very early stage of a 
heavy-ion collision. 

We obtain a rather interesting picture of the very early 
stage of ultrarelativistic heavy-ion collisions. Magnetic 
Wilson loops of area > 2/Ql exhibit area law behavior 
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FIG. 4: Area enclosed by a quantum mechanical path shifted 
by about one de Broglie wavelength from the extremal classi- 
cal path. 



which implies uncorrelated magnetic vortex-like flux be- 
yond the scale R^tx ^ 0.8/ Qs- This indicates that many 
more multipole moments rather than just the dipole are 
present. The vortex structure of the longitudinal mag- 
netic field modifies propagation of particle-like modes 
with de Broglie wavelength somewhat larger than Qs. 
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